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Abstract 

The error between appropriately smooth functions and their radial basis function 
interpolants, as the interpolation points fill out a bounded domain in H'', is a well 
studied artifact. In all of these cases, the analysis takes place in a natural function 
space dictated by the choice of radial basis function — the native space. The native 
space contains functions possessing a certain amount of smoothness. This paper es- 
tablishes error estimates when the function being interpolated is conspicuously rough. 
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1 Introduction 



In this paper we are interested in interpolation of a finite scattered data set A C M*^ by 
translates of a single basis function. Of the differing set ups to this problem, the one 
preferred in this paper is the following variational formulation. Firstly, we require a space 
of continuous functions Z which carries a seminorm. The minimal norm interpolant to 
/ : — > H on ^ from Z is the function Sf € Z which agrees with / on ^ and has 
smallest seminorm amongst all other interpolants to / on ^ from Z. The particular space 
we shall be concerned with is 

^'-(]R'^):=|/G^': D^f G Li,ioe(]R''), J^^ w (x) \ (D^f )ix)\^ dx < oo, |a|=m|, 



which carries the seminorm 

|/U:=(^ c^j ^w{x)\{D^f){x)\^dx\"\ /G2'"(]R'^). 



The constants are chosen so that 'Yl,\a\=m^ax'^'^ = for all x G JM^. The notation 

=5^' is used to denote the usual Schwartz space of distributions. The space Z^{M!^) is 
christened the native space. The weight function function w : M"^ ^ H is initially chosen 
to satisfy 

(WO) w G C7(]R'^\0); 
(Wl) w{x) > if X / 0; 
(W2) l/w G Li,ioe(Il''); 

(W3) there is a positive /i G IR such that {w{x))~^ = 0{\x\~'^^^) as |x| oo. 

A consequence of (W0)-(W3) is that Z™(IR'^) is complete with respect to | • U, and if 
m + — d/2 > then Z™'(R'^) is embedded in the continuous functions (see [6j). As the 
title of this work suggests, we expect this set up to admit minimal norm interpolants of 
the form 

(5/)(x) := ^6,V(x-a), forxGlR^ (1.1) 

for an appropriate basis function ip. We are not disappointed, but for brevity we omit the 
details which are well presented in The coefficients ba in are determined by the 
interpolation equations {Sf){a) = f{a), a ^ A. In some situations it may be necessary 
to append a polynomial p onto p.l|) and take up the ensuing extra degrees of freedom by 
satisfying the side conditions: 

^ baq{a) = 0, 

a&A 

whenever g is a polynomial of the same degree (or less) as p. The archetypal scenario 
the author has in mind is w{x) = \x\'^^ for x G M'^, where ^ < d/2. This leads to 
minimal norm interpolants of the form (jl.ip modulo a polynomial of degree m. Here, the 
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radial basis function is il) : x t-^ |^|2m+2^~d j^| •£ 2m + 2fi — d is an even integer or 
■0 : X I— > |xp™+2/^-'^ otherwise. 

It is of central importance to understand the behaviour of the error between a function 
/ : $7 — > ]R and its interpolant as the set A becomes dense in a bounded domain J7. The 
measure of density we employ is the fill-distance h := sup^g^ miuaeyilx — a\. One finds 
that there is a positive constant 7(717.) , independent of h, such that for all / S Z"^{1R'^), 

\\f-Sf\\L,in) = 0{h^'~'^^), ash^O. 

It is natural to ask: what happens if the function being approximated does not lie in 
2™(]R'^)? It may well be that / lies in 2:^(11'^), where k < m and k + ^ - d/2 > 0. The 
condition k + fj, — d/2 > ensures that /(a) exists for each a G so Sf certainly exists. 
It is tempting to conjecture that the new error estimate should be 

\\f-Sfh,in) = Oih^^''^), ash^O. 

We are conjecturing the same approximation order as if we had instead approximated 
/ with the minimal norm interpolant to / on ^ from ^'^(IR'^). This is precisely what 
happens in the case w = 1, which was considered by Brownlee &: Light in [_2j. In this 
work, with the aid of a recent result from pj (Lemma [23]), we employ the technique used 
by Brownlee & Light to extend their work to more general weight functions. Theorem 13.51 
is the definitive result we obtain. The interested reader may enjoy consulting the related 
papers [3 El EJ HO] . 

To close this section we introduce some notation that will be employed throughout the 
paper. A domain is understood to be a connected open set. The support of a function 
^ : R"^ ^ M, denoted by supp {(f>), is defined to be the closure of the set {x G M'^ : (f){x) 7^ 
0}. We make much use of the linear space 11^(11'^) which consists of all polynomials of 
degree at most m in d variables. We fix i as the dimension of this space. Finally, when we 
write / we mean the Fourier transform of /. The context will clarify whether the Fourier 
transform is the natural one on Li(lR'^), f{x) := (27r)^'^/^ J^^d /(t)e~'^* dt, or one of its 
several extensions to L2(1R'^) or J/" . 

2 Extension theorems 

In this section we gather a number of useful results, chiefly about the sorts of extensions 
which can be carried out on our native spaces. This will first require us to establish 
the notion of local native spaces. To do this, we rewrite the seminorm \f\m in its direct 
form — that is, without the Fourier transform of / appearing explicitly. Let us demand 
that w satisfies the following additional axioms: 

(W4) w{y) = w{-y) for all y G R''; 

(W5) ty(0) = and w{x)<Q for almost all x G R''; 

(W6) w isa measurable function and for any neighbourhood of the origin, w G L^(1R'^\ 
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(W7) \w {y)\ = 0{\y\^) as y ^ , where A + d + 2 > 0. 

Armed with axioms (Wl) and (W4)-W7) it follows from [3] that 

\f\l = -lT.cJ' [ wix-y)\iDy){x)-iDy){y)\'dxdy, / € Z-(]R'^). 

(2.1) 

The notation C^(IR'^) is used for the space of compactly supported m-times continuously 
differentiable functions on M'^. Now, let us define the following space for a domain $7 C K*^, 



where 



\f\m,n:=(-l Yl / i^ix-y)\iDV)ix)-iDy)iy)\''dxdyy\ / G 

A norm is placed on via 

||/IU,Q := {WfWw^^in) + \f\l,ny, f e ^™(^^)- 

The notation denotes the completion of X"^{0,) with respect to || • \\m.n, while 

y^^{Q) denotes the completion of with respect to | • \m,n- It is these spaces that 

we call the local native spaces. 

We are nearly ready to state our first extension theorem, but first it is necessary to take 
on board four additional axioms and introduce an important type of bounded domain: 

Definition 2.1. Let and domains in TR'^, and ^ a bijection from to ^2- 

We say that <I> is m-smooth if, writing (^{x) = {(pi^xi, ... ,Xd), ■■■ jCpdixi, ■■■ ,Xci)) and 
^^^{x) = ^'(a;) = {ipi{xi, . . . , Xd), . . . , ipdixi, ■ ■ ■ , Xd)), then the functions (/>!,..., 0^ belong 
to C"^(r2i) and V'l, • • • , V'd belong to C"(r22)- Let ^ be a bijection from M!^ to TR^. We 
say $ is locally m-smooth if ^ is m-smooth on every bounded domain in M!^. 

(W8) for every locally (m + l)-smooth map (p on R'', and every bounded subset O of H , 
there is a Ci > such that w {cj){x) — 4>{y)) < Ciw [x — y), for all x,y £ Q; 

(W9) there exists a constant C2 > such that if x = {x' , Xd) G M'^ and y = {x' , yd) G K*^ 
with \xd\ > \yd\, then id (x) < C2W (y). 

(WIO) J^tv < whenever A has positive measure; 

(Wll) w{y)=w i-y) for ah y G W^. 

Definition 2.2. Let B = {(2/1,2/2, • • • , yd) G R'^ : \yj\ < 1, 1 < j < d}, and set = 

{y e B : y = {y', yd) and yd > 0} and Bq = {y e B : y = (y', y„) and yn = 0}. A bounded 
convex domain in M!^ with boundary will be called a V-domain if the following all 
hold: 
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(Al) there exist open sets Gi, . . . ,Gm C K*^ such that C UjLi Gj] 

(A2) there exist locally (m+lj-smooth maps cj)j : H'^ H"^ such that (t)j{B) = Gj, 
(j)jiB+) = Gjnn and (j)j{Bo) = Gj n dQ, j = I, . . . , N; 

(A3) let 0,s be the set of all points in 0, whose distance from is less than 5. Then 
for some 6 > 0, 

5^5 C [J 0j N (yi, 2/2, . . . , yd) G H'^ : \yj\ < ^<j<d 
j=i 

The definition of a V-domain is taken from a paper by Light & Vail [3] in which 
extension theorems for our local native spaces are considered. 

Theorem 2.3 (Light k Vail [5j). Let n C H'^ be a V-domain. Let w : R'' ^ IR satisfy 
(W6)-(W11). Then there exists a continuous linear operator L : — > ^'"(IR'^) such 
that for all f G X'^{VL), 

1. Lf = f on ri; 

2. supp {Lf) is compact and independent of f; 

3. \\Lf W^-^d < i^||/||m,Q, for some positive constant K = K{^) independent of f. 

A feature of the construction of the extension operator in Theorem 12.31 is that Lf can 
be chosen to be supported on any compact subset of IR'^ containing 17. For details of the 
construction, the reader should consult [5j. Also at our disposal is a seminorm version of 
Theorem 12.31 

Theorem 2.4 (Light k Vail [S]). Let n C H'^ be a V-domain. Let w : R'^ ^ IR satisfy 
(W6)-(W11). Given f G 3^™(J7), there exists a function f^ G 3^™(1R'^) such that: 

1. f'' = f on n; 

l^^lmR'' — G\f\m,n, for some positive constant G = C{Q) independent of f . 
It is convenient for us to be able to work with a norm on ^™(il) that is equivalent to 

II ■ ||m,f2 • 

Lemma 2.5 (Brownlee k Levesley [Tj). Let Cl C R"^ be a V-domain. Let : R'^ ^ R 
satisfy (W0)-(W12) and let m + ^ — d/2 > 0. Let bi, . . . ,bi £ Q be unisolvent with respect 
to nm(R'^). Define a norm on A'™' ($7) via 

^ 1=1 

There are positive constants Ki and K2 such that for all f G X'^{Vl), i^i||/||m,c ^ ll/llo < 

K2\\f\\m,n. 
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The behaviour of the constant K{i}) in the statement of Theorem l2.3l can be understood 
for simple choices of 0. To reahse this, we require that the weight function satisfies one 
further and final axiom: 

(W12) there exists Ci,C2 > such that Cih^w (x) < w {hx) < C2h^ w (x), for all /i > 0, 
X G R'^. 

Now, an elementary change of variables gives us: 

Lemma 2.6. Let J7 6e a measurable subset of TR'^. Let w : IR'^ — > IR 6e a measurable 
function that is nonpositive almost everywhere and satisfies (Wll). Define the mapping 
a : — > H"^ by ct(x) = a + h{x — t), where h > 0, and a, t, x G IR'^. Then there exists a 
constant Ki, K2 > 0, independent of fi, such that for all f E y"^(a{0,)), 

^ \focr\m,n ^ 

^1 ^ h-^f^-VU,,, - 

We are now ready to state the key result of this section, but before doing this let us 
make a simple observation. Look at the unisolvent points 61, . . . , 6^ in the statement of 
Lemma 12.51 Since X"^(i}) can be embedded in C(0), it makes sense to talk about the 
interpolation operator P : — > UmiM!^) based on these points. 

Lemma 2.7. Let w : IR^ IR satisfy (W0)-(W12). Let B be any ball of radius h and 
centre a G IR'^, and let f G X'^{B). Whenever 61, . . . , 6^ G IR'^ are unisolvent with respect 
to nm(IR'^) let Ph : C(1R'^) — 11^(11'^) be the Lagrange interpolation operator on bi, . . . ,b£. 
Then there exists c = (ci, . . . , q) G i?^ and g G X'^{1R'^) such that 

L g{x) = {f-PJ){x)forallxeB; 

2. g{x) = for all |x - a| > 2h; 

3. there exists a C > 0, independent of f and B, such that \g\^jid < C\f\m,B- 
Furthermore, ci, . . . , q can be arranged so that ci = a. 

Proof. Let Bi be the unit ball in TR'^ and let B2 = 2Bi. Let 61, . . . , 6^ G -Bi be unisolvent 
with respect to nm(lR'^)- Define a{x) = — a) for all x G IR*^. Set q = a~^{bi) 

for i = 1,...,^ so that ci,...,C£ G B are unisolvent with respect to Ilm{^'^)- Take 
/ G Then (/ - PJ) o a'^ G X"^{Bi). Set F = {f - PJ) o a-\ Let F^i be 

constructed as an extension to F on Bi. By Theorem 12.31 and the remark following it, we 
can assume F^^ is supported on B2. Define g = F^^ o a G X^{M!^). Let x ^ B. Since 
(j{B) = Bi there is a y G i?i such that x = a~^{y). Then, 

g{x) = (F^^ o a){x) = F^^{y) = ((/ - PJ) o a-'){y) = (/ - PJ){x). 

Also, for X G K'^ with |rE — a| > 2/i, we have \cr{x)\ > 2. Since F^^ is supported on B2, 
g{x) = for la; — a| > 2h. Hence, g satisfies properties 1 and 2. By Theorem 12.31 there is 
a Ki, independent of / and B, such that 

\\F^'\\m,B2 < II^^ML,Rd < Ki\\F\\^,B,. 
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We have seen in Lemma 12.51 that if we endow X"^{Bi) and X™'(B2) with the norms 




then 11 ■ 11^. and || • \\m,Bi are equivalent for i = 1,2. Thus, there are constants K2 and K3, 
independent of / and B, such that 

WP'^'We, < K2\\F'''\\m,B, < KiK2\\F\\m,B, < KiK2K^\\F\\b,. 

Set Ci = KiK2K3. Since F^^{bi) = F{bi) = (/ - PJ){a-\bi)) = (/ - Pcf){ci) = for 
i = !,...,£, it follows that \F^'\m,B2 < Ci\F\m,B^. Thus, \g o a'^l^^^, < Ci\{f - PJ) o 
o'~^\m,Bi- Now, Lemma 12.61 can be employed twice to provide us with constants C2 and 
C3 > 0, independent of / and B, such that 

< CiC2h''+^'^-^\{f - PJ) O fT-V,Bi < C1C2C3I/ - Pcf\m,B- 

Finally, we observe that \ f — Pcf\m,B = \f\m,B to complete the first part of the proof. The 
remaining part follows by selecting 61 = and choosing b2, ■ ■ ■ ,bi accordingly in the above 
construction. □ 



3 Error estimates 

In this section we establish the error estimate conjectured in the introduction. We begin 
with a function / in Z''(R'^). We want to estimate 

\\f - Smfh^iU), (3.1) 

where Sm is the minimal norm interpolation operator from Z"^(R'^) on A and m > k. 
The essence of the proof is as follows. Firstly, by adjusting /, we obtain a function /, 
still in Z^'QR'^), with seminorm in Z^(JR'^) not too far from that of /. We then smooth 
/ by convolving it with a function (f) G C^(]R'^). The key feature of the adjustment of 
f to F := 4> * f is that F{a) = f{a) for every point a € A (Theorem 13. 4p . This enables 
us to replace Smf with SmF in (13. ip . Furthermore, it follows that F G Z'^{1R'^) so we 
can employ an existing L2-error estimate to -F — SmF. The remaining part of the error, 
f — F, is easily dealt with as it vanishes on A. Finally, Lemma 13.11 takes us back to an 
error estimate in 2*^(11'^). 

Lemma 3.1. Let w : Wif^ ^ TR satisfy (WO) and (Wl). Let k < m and (p G C^(Il'^). 

For each h > let (phix) = h~'^(j){x/h) for x G H'^. Then there exists a constant C > 0, 
independent of h, such that for all f G Z''{H'^), \(j)h * flmn'' — ^^''~^\f\kn'^- 

Proof. The case w = 1 is established in [2|. The proof for this more general set up does 
not differ substantially so is omitted. □ 
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Lemma 3.2 (Brownlee & Light |2j). Suppose 4> S C^{1R'^) is supported on the unit ball 
and satisfies 

/ (f>{x) dx = 1 and / (f){x)x°' dx = 0, for all < |a[ < k. 

For each e > and x G TR'^, let (f)e{x) = e~'^4'{x/£). Let B he any hall of radius h and 
centre a G H'^. For a fixed p G Yli^iJM^) let f he a mapping from H"^ to M such that 
f{x) = p{x) for all X £ B. Then {(p^ * = p{a) for all e < h. 

Definition 3.3. Let ri he an open, bounded subset o/M'^. Let A be a set of points in 
The quantity h := sup^g^ infag^|x — a| is called the fill-distance of A in 17. The separation 

of A is given by the quantity q := miUa^beA ^'~2~^ ■ The quantity h/q will he called the 
mesh-ratio of A. 

Theorem 3.4. Let w : ^ U satisfy (W0)-(W12). Let k + fi - d/2 > and m> k. 

Let A be a finite subset o/H^ of separation q > 0. Then for all f G X^(iR'^) there exists 
anF e X"^(R'^) such that 

1. F{a) = f{a) for all a G A; 

2. there exists a C > 0, independent of f and q, with jj^d < C|/|^ and \F\^ < 

Proof Take / G For each a G ^ let C R'^ denote the ball of radius 6 = q/4 

centred at a. For each Ba let ga be constructed in accordance with Lemma 12.71 That is, 
for each a £ A take c' = (c2, . . . , q) G B^^^ and ga G X''{1R'^) such that 

1. a, C2, . . . , Q are unisolvent with respect to 11^(11'^) 

2. ga{x) = if - P(a,c')f){x) for ah x G Ba] 

3. Pi^a,c')f G nfc(Il'^) and (P(a,c')/)(a) = /(«); 

4. ga{x) = for all |a; - a| > 25; 

5. there exists a Ci > 0, independent of / and Ba, such that balfc^d < C'i!/|fe,_Ba- 

Note that if a 7^ 6, then supp(5(a) does not intersect supp((7fe), because if x G supp(5(a) 
then 

|x — 6| > |6 — a| — |x — a| > 2(/ — 2(5 = 6(5. 
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Let U = UfcG^ supp {gt), then writing R"^ = (R'^ \U)UU we obtain 



2 I 



w (x — y) 



2 ^ 

\a\=k 
+ 2 



w{x-y) 



w {x — y) 



dxdy 



^((D"5a)(x)-(Z)V)(y)) 



dxdy 



^((Z)"5a)(x)-P"5a)(y)) 



+ / I w{x-y) 
Ju Ju 



dxdy 



^iiD-ga)ix)-iD-gM) 
aeA 



dxdyj . 



(3.2) 



We shall now consider each of the double integrals in (|3.2p separately. Firstly, the integral 
over (R'^ \U) x (K*^ \ U) is zero because Yla&A9<^ supported on U. Next, using the 
observation above regarding the support oi ga, a G A, it follows that 

2 



4E 



w {x — y) 



Y^{{D^ga){x) - {D^ga){y)) 



dxdy 



h&A H=k -J^^KU-Jsuppia,) 



w{x-y) 



a£A 



dxdy 



e4e 



b&A \a\=k JTR^XU Jsupp(g,) 



wix-y)\iD'^gb)ix)\^dxdy 



= T.-l E 

6G^ \a\=k 

< Elff"! 



w{x-y)\iD''gb)ix) - iD''gb)iy)\^ dxdy 



2 

kM''' 



(3.3) 



b€A 



Before calculating the final integral let us examine the following expression for b € A and 
a G TL^ with lal = m, 



E 



cg_4 -'supp (gc) Jsupp (gb) 

cGA •^'^UPP (Sc) -'supp (gt) 
c^b 

^2E 

CG.A -^supp (gc) ./supp (g;,) 



(x - y) I (D'^gt,) (x) - {D^g,) {y) ^ dxdy 

w{x-y)\{D^g,){x)- {D^g,) (y) + {D^g,) {x) - {D^ g,){y)\'' dxdy 

w{x- y)\{D''g,){x) - {D^gb){y)? dxdy 

2J2[ [ w{x-y)\{D''g,)ix)-iD^g,){y)\^dxdy 

ceA '^"PP ^^"^ '^"PP ^Sb) 

c^b 
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< 2 



w{x- y)\{D''g^){x) - {D'^9,){yt dxdy 



w{x- y)\{D^g,){x) - {D^9c){y)? dxdy. (3.4) 



Finally, using the observation regarding the support of (^a, a ^ A, once again and ()3.4|) it 
follows that 



4e 



\a\=k 



w{x-y) 



u JU 



Y,{{D^9a){x)-{D^ga){y)) 



dxdy 



beAceA \a\=k J ^^PP idc) J supp (g^) 



w {x — y) 



Y,{{D^9a){x)-{D-ga){y)) 
aeA 



dxdy 



bdAceA \a\=k -J ^^PP {9c) J suppigt) 



iv{x-y)\ {D^gb){x) - {D^g,) (y) ^ dxdy 



+ee4e , 

be^ c£A 

^E 



u5 (x - 2/) I (x) - (y) I' dxdy 



w}(x - y)|p"5fe)(a;) - P"<7c)(y)P dxdy 



2 



< 



Substituting and into ([321) we find 



(3.5) 



Y.9a 

aeA 



^ - '^X]l5'«lfc,Rd- 



Hence, applying Condition 5 to the above inequality we have 

Y.9a' ^<7CfJ2\f\ls^<7C!\f\ 
aeA aeA 



2 

fc.R'*' 



Now set H = f — YlaeAd"-- then follows from Condition 1 that H{x) = {P(a,c')f ){x) for 

'0 



all X e Ba, and from Condition 3 that H{a) = f{a) for all a G ^. Let e C^(1R'^) be 



supported on the unit ball and enjoy the properties 



/ 



(x) dx = 1 and 



(x)x" dx = 0, for all < lal < k. 



R" 



Now set F = (ps * H . Using Lemma |3.H there is a constant C2 > 0, independent of q and 
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/, such that 



<2C2<52(fc-™)(|/|2^,+ 



<2C2{l + lCl)6'- i.i^jj^,. 

Similarly, there is a constant C3 > 0, independent of q and /, such that 

2 



<C3 



9a 



< 2^3(1 + 7C7f)|/| 



2 



Thus \F\^-^d < Cq and \F\j^^d < C\f\f.^d for some appropriate constant 



C > 0. Since F = (ps * H and i^l^^E 11^(11 ) for each a S it follows from Lemma 
that F{a) = H{a) = f{a) for all a e A. □ 

Theorem 3.5. Let C H'^ he a V-domain and let w : IR'^ ^ IR, 6e a measurable function 
satisfying (W0)-(W12). Let k + fi - d/2 > and m > k. For each h > 0, let Ah be 
a finite, Ilm(R'^)-unisolvent subset of fl with fill- distance h. Assume also that there is a 
quantity p > such that the mesh-ratio of each Ah is bounded by p for all h > 0. For each 
mapping f : Ah H, let S^f be the minimal norm interpolant to f on Ah from Z'^(R'^). 
Then there exists a constant C > 0, independent of h, such that for all f S y^{^), 

11/ - Sif h,(^n) < Ch''^/'-''/V\k,n, ash^O. 

Proof. Take / G X''(R'^). Construct F in accordance with Theorem [331 and set G = f — F. 
Then F(a) = f{a) and G{a) = for all a G Ah- Furthermore, there is a constant Ci > 0, 
independent of / and /i, such that 



\F\ 



h\k- 



Thus Sti 

Hence, 



StF and S'^G 



\f\k,u'^ \G\k,u^ < \f\k,u^ + \F\k^ud < (1 + Gi)\f\^^^d. (3.6) 
0, where we have adopted the obvious notation for S^. 



11/ - 'S'm/llL2(n) -\\{F + G) - S^F\\l^(p^^ < \\F - S'^F\\i^^q) + ||G - S^GUl^^^^). 

Now, employing the error estimate in \T], there are positive constants C2 > and C3 > 0, 
independent of h and /, such that 

11/ - Stf h,in) < C2h"'~'^'-''/^\FUn + G3h''-'/'-'''/^\G\k,n, as h ^ 0. 
Finally, using the bounds in ()3.6p we have 

11/ - Sifh,in) < C^h'^-'^'-'^/VlkM''^ ash^O, (3.7) 



for some appropriate C4 > 0. In particular, ((TTD holds for all / G X^(IR'^). As y''(R'^) 
is a dense linear subspace of X''{TR'^) then (13.70 extends to hold for all / G y^{lR'^) using 
a standard normed space argument [31 Page 180]. To complete the proof we now let 
/ G y^{Q) and define in accordance with Theorem 12.41 It follows that there is a 
C5 > such that 

11/ - Sif h,in) < C,h'-''^^-''/^\f\n<i < C,C5h''-^/'-''/^\f\k,n, as h ^ 0. □ 
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